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The development of crops which are well suited to growth 
under future environmental conditions such as higher at- 
mospheric CO2 concentrations ([CO2]) is essential to meet- 
ing the challenge of ensuring food security in the face of the 
growing human population and changing climate. A high- 
yielding indica rice variety {Oryza satlva L. cv. Takanari) has 
been recently identified as a potential candidate for such 
breeding, due to its high productivity in present [CO2]. To 
test if it could further increase its productivity under ele- 
vated [CO2] (eC02), Takanari was grown in the paddy field 
under season-long free-air CO2 enrichment (FACE, approxi- 
mately 200 |imol mol~^ above ambient [CO2]) and its leaf 
physiology was compared with the representative japonica 
variety 'Koshihikari'. Takanari showed consistently higher 
midday photosynthesis and stomatal conductance than 
Koshihikari under both ambient and FACE growth condi- 
tions over 2 years. Maximum ribulose-l,5-bisphosphate 
carboxylation and electron transport rates were higher for 
Takanari at the mid-grain filling stage in both years. 
Mesophyll conductance was higher in Takanari than in 
Koshihikari at the late grain-filling stage. In contrast to 
Koshihikari, Takanari grown under FACE conditions 
showed no decrease in total leaf nitrogen on an area basis 
relative to ambient-grown plants. Chi content was higher 
in Takanari than in Koshihikari at the same leaf nitro- 
gen level. These results indicate that Takanari maintains 
its superiority over Koshihikari in regards to its leaf-level 
productivity when grown in elevated [CO2] and it may be 
a valuable resource for rice breeding programs which seek 
to increase crop productivity under current and future 
[CO2]. 



Keywords: Ecophysiology 
Photosynthesis. 



• FACE • Global change • 



Abbreviations: A, net photosynthetic carbon assimilation 
rate; C^, atmospheric CO2 mole fraction; CO2 mole frac- 
tion at the site of carboxylation; CO2 mole fraction in the 
intercellular airspace; eC02, elevated [CO2]; FACE, free-air 
CO2 enrichment; GF, grain filling; g^, mesophyll conductance 
to CO2; gs> stomatal conductance to water vapor; J^ax' the 
maximum electron transport rate supporting the regener- 
ation of ribulose-1,5-bisphosphate; MCL, mesophyll conduct- 
ance limitation to photosynthesis; N, nitrogen; Rubisco, 
ribulose-1,5-bisphosphate carboxylase oxygenase; SL, stomatal 
limitation to photosynthesis; V^, max/ the maximum carboxyl- 
ation rate of Rubisco. 



Introduction 

As the world population continues to increase in the coming 
decades, the demand for food is expected to increase substan- 
tially (Alexandratos and Bruinsma 2012). However, the progres- 
sion of climate change and an increasing frequency of extreme 
weather events in various parts of the world are expected to 
have negative effects on food crop yields in the future (Hatfield 
et al. 2011, Yadav et al. 2011). To ensure food security going 
forward, it is essential that we gain an understanding of how 
crops respond to the changing environment and use that in- 
formation to develop new varieties of staple food crops which 
can tolerate future conditions or even take advantage of them 
to become more productive than current cultivars (Ainsworth 
et al. 2008, Ziska et al. 2012, Tausz et al. 2013). 

The rise of atmospheric concentrations of CO2 ([CO2]) since 
the Industrial Revolution is one of the most clearly documented 
aspects of global change, and the global mean [CO2] is projected 
to continue to rise over the course of this century (Meehl et al. 
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2007). The primary effect of elevated [CO2] (eC02) on the physi- 
ology of C3 crops is an increase in the rate of photosynthesis, 
although this does not always translate into an equivalent 
increase in biomass accumulation or yield (Leakey et al. 2009). 

Rice (Oryza satwa L.) is the staple food for over half of the 
world's population, and by 2050 the demand for rice is ex- 
pected to increase by nearly 30% over 2005-2007 production 
levels (Alexandratos and Bruinsma 2012). The genetic diversity 
of rice is high, and improving its yield potential via selective 
breeding of key cultivars will be a large step towards achieving 
higher productivity in the future (Khush 2005). 

A high-yielding \nd\ca variety of rice called 'Takanari' was 
developed in Japan in the 1980s as the offspring of two high- 
yielding cultivars from Korea ('Milyang 25' and 'Milyang 42') 
(Imbe et al. 2004, Takai et al. 2012). Compared with the average 
japonica rice variety, Takanari exhibits some remarkable physio- 
logical characteristics, such as a very high stomatal and 
hydraulic conductance, as well as a large below-ground rooting 
system which enables it to accumulate a large amount of ni- 
trogen (N) (Taylaran et al. 2009, Taylaran et al. 2011). These 
physiological traits result in high source capacity (i.e. leaf-level 
photosynthesis and carbon assimilation) (Taylaran et al. 2011) 
and, coupled with a high sink capacity (i.e. large panicle size and 
excellent grain-filling efficiency) (Nagata et al. 2001, Hirasawa 
et al. 2010), are responsible for its high biomass and yield when 
grown under present-day field conditions. Unfortunately, the 
quality and edibility of its grain are low (Imbe et al. 2004), but it 
remains a potential germplasm resource for breeding 
higher productivity into the next generation of modern rice 
varieties. 

Consequently, to examine its response under future atmos- 
pheric conditions in the field, Takanari has been grown at two 
free-air CO2 enrichment (FACE) experimental sites in Japan 
under a season-long open-air fumigation of elevated [CO2] 
(+200 |imol moP^ above ambient [CO2]) (Hasegawa et al. 
2013). FACE technology involves the computer-controlled 
release of CO2 from an array of pipes or blowers laid out in 
the field; the gas is then carried across the treatment area via 
natural wind and diffusion (Hendrey and Miglietta 2006). It has 
been used for >20 years to investigate the response of a variety 
of natural and managed ecosystems to eC02 (reviewed by 
Ainsworth and Long 2005, Nosberger et al. 2006). 

At the Shizukuishi FACE site in northeastern Japan during 
the 2008 growing season, Takanari showed a 16% yield en- 
hancement in eC02, compared with 10% for 'Koshihikari' (cur- 
rently considered the representative japonica variety, and 
grown widely across Japan) (Hasegawa et al. 2013). In 2010, at 
the Tsukubamirai FACE site in central Japan, Takanari showed a 
21% yield enhancement in eC02 compared with 16% for 
Koshihikari (Hasegawa et al. 2013). The brown rice yield 
(g m~^) of Takanari was 26% and 32% higher than that of 
Koshihikari under ambient and eC02 conditions, respectively, 
indicating that the large sink capacity of Takanari (specifically, 
high grain number per panicle and high panicle number per 
hill) was effective in enabling a strong yield enhancement 



response in eC02 (Hasegawa et al. 2013). However, the response 
of the source capacity (i.e. photosynthesis) and the leaf-level 
physiology of Takanari under eC02 have not yet been 
examined. 

The objective of this study was to characterize and compare 
the leaf-level physiology of the rice cultivars Takanari and 
Koshihikari under a season-long, free-air elevated [CO2] treat- 
ment. We sought to answer the question of whether the source 
capacity of Takanari during grain filling is greater than that of 
Koshihikari in eC02, by examining the photosynthesis, the sto- 
matal and mesophyll conductance, leaf N, protein content and 
pigment content of the two cultivars from the booting stage to 
the late grain-filling (GF) stage over parts of two growing seasons. 




Photosynthesis and stomatal conductance 

The photosynthetic rate (A) of Takanari flag leaves was signifi- 
cantly higher than that of Koshihikari in both the ambient and 
eC02 treatments across all measurement time points in 2012 
and 2013, except for the ambient-grown plants at the booting 
stage (1.5 weeks before heading) in 2012 (Fig. 1A; Table 1). The 
highest photosynthetic rate for Takanari flag leaves in the FACE 
treatment (at 590|imol moP^) was 44.2 |imol CO2 m~^ s~\ 
observed at the early GF stage (1 week after heading) in 2012. In 
comparison, the maximum observed midday photosynthesis of 
Koshihikari in FACE was 36.4 |imol CO2 m~^ s~^ at the mid GF 
stage (2 weeks after heading) in 2012. The relative advantage of 
Takanari over Koshihikari at each growth stage was generally 
similar between the ambient and FACE treatments. For ex- 
ample, during the booting stage in 2012, A of Takanari was 
34.0% and 34.2% higher than that of Koshihikari in the ambient 
and FACE treatments, respectively. This advantage in Takanari 
decreased over time to 14.0% and 14.4% at the late GF stage 
(4 weeks after heading) in 2012, resulting in a significant cultivar 
by stage effect for A in 2012 (P = 0.046, Table 1). In 2013, A in 
Takanari was consistently higher than in Koshihikari from boot- 
ing to mid GF (2 weeks after heading) in both ambient and 
FACE treatments, but the highest advantage of Takanari was 
observed at early GF (1 week after heading) (45.7% and 42.7%, 
respectively). 

The overall effect of FACE on photosynthesis was significant 
in both years (Fig. 1A; Table 1). The relative enhancement of A 
in plants grown in FACE over ambient conditions was compar- 
able in the two cultivars within each growth stage. In 2012, the 
enhancement of A by FACE at the booting stage in Koshihikari 
and Takanari was 35.6% and 35.8%, respectively; this decreased 
to 19.3% and 19.7% at the late GF stage. During the 2013 field 
season, the enhancement of A by FACE ranged from 17.1% to 
35.0% in Koshihikari and from 19.5% to 30.9% in Takanari. 

Takanari consistently showed a significantly higher stomatal 
conductance (gs) than that of Koshihikari, regardless of the 
growth stage or growth [CO2] in both years, with the exception 
of the booting stage in 2012 (Fig. IB; Table 1). The difference in 
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Fig. 1 Midday photosynthetic gas exchange measurements taken at the respective growth [CO2] for ambient (plain bars) and FACE (cross- 
hatched bars) treatments (390 and 590 |imol mor\ respectively). Koshihikari (KH) is shown in the white bars and Takanari (TN) in the gray bars. 
Data are from the 2012 and 2013 field seasons. The measurements were taken at the following growth stages: (1) booting (approximately 
1.5 weeks before heading); (2) full heading (3 d after heading 2013 only); (3) early grain filling (GF) stage (1 week after heading); (4) mid GF stage 
(2 weeks after heading); and (5) late GF stage (4 weeks after heading 2012 only). All measurements were taken on flag leaves. Error bars are SEMs. 
Different letters above each bar indicate significant differences between means at a = 0.05 within each measurement time point. n = 4. (A) Net 
photosynthetic carbon assimilation rate (A). (B) Stomatal conductance to water vapor (gj. 



gs between cultivars varied between measurement time points, 
butgs reached as high as approximately 100% and 195% higher 
in Takanari than in Koshihikari at early GF in 2012 and 2013, 
respectively. There was a significant overall effect of FACE on gs 
over the entire experiment in both years (P = 0.019, 0.0003; 
Table 1), but this FACE-induced decrease in gs was only signifi- 
cant at the level of individual treatment mean comparisons in 
Takanari at the late GF in 2012 and at booting and full heading 
(3d after heading) in 2013 (Fig. IB). 

Biochemical capacity for photosynthesis 

Takanari showed a significantly higher biochemical capacity for 
photosynthesis than Koshihikari, as indicated by higher values 
^cmax [the maximum carboxylation rate of ribulose-1, 
5-bisphosphate carboxylase oxygenase (Rubisco)] and j^ax 



(the maximum electron transport rate supporting the regener- 
ation of ribulose-1,5-bisphosphate) (Fig. 2A, B; Table 1). V^max 
and Jmax were both higher in Takanari than in Koshihikari at the 
full heading stage in 2013 and mid GF in 2012 and 2013. 
However, at late GF in 2012, there was no significant difference 
•ri Vcmax and Jmax between either cultivars or CO2 treatment. 

CO2 diffusion-related limitations to 
photosynthesis 

There was a significant overall effect of cultivar, growth [CO2] and 
stage on mesophyll conductance (gm) in both 2012 and 2013 
(Table 1). During the 2012 growing season, Takanari had a 
higher g^ than Koshihikari in the FACE treatment at mid GF 
and in both ambient and FACE treatments at late GF (Fig. 2C). 
In 2013, gm was identical between cultivars at full heading but was 
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Table 1 Three-way analysis of variance of the midday net photosynthetic carbon assimilation rate (A), stomatal conductance to water 
vapor (gs), maximum carboxylation rate of Rubisco iVc,max)> maximum electron transport rate (Jmax)> mesophyll conductance to CO2 (gmX 
stomatal limitation to photosynthesis (SL) at each treatment's respective growth [CO2] and mesophyll conductance limitation to photo- 
synthesis (MCL) at each treatment's respective growth [CO2] from Figs. 1 and 2 



Effect 


A 


gs 


^c,max 


J max 




SL 


MCL 


2012 growing season 
















Cultivar 


<0.0001 


<0.0001 


0.0033 


<0.0001 


0.0002 


0.0003 


0.0005 


Growth CO2 


<0.0001 


0.0191 


NS 


0.0134 


0.0018 


<0.0001 


<0.0001 


Stage 


<0.0001 


<0.0001 


<0.0001 


0.0321 


<0.0001 


NS 


<0.0001 


Cultivar x CO2 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


Cultivar x Stage 


0.0457 


NS 


0.0207 


NS 


NS 


NS 


0.0084 


CO2 X Stage 


0.0498 


NS 


NS 


NS 


NS 


NS 


NS 


Cultivar x CO2 x Stage 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


2013 growing season 
















Cultivar 


<0.0001 


<0.0001 


<0.0001 


0.0117 


0.0005 


0.0005 


0.062 


Growth CO2 


<0.0001 


0.0003 


NS 


0.0021 


0.0023 


0.0002 


<0.0001 


Stage 


0.0039 


<0.0001 


0.0149 


<0.0001 


0.0003 


<0.0001 


NS 


Cultivar x CO2 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


Cultivar x Stage 


NS 


<0.0001 


NS 


0.087 


0.0025 


NS 


NS 


CO2 X Stage 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


Cultivar x CO2 x Stage 


NS 


NS 


NS 


NS 


NS 


NS 


0.082 



Data are from the 2012 and 2013 growing seasons. 

P-values <0.05 are in bold; P-values >0.05 but <0.10 are in italics; P-values >0.10 are indicated by NS. 



Table 2 Three-way analysis of variance of the time course data for Rubisco content, leaf soluble protein, total leaf nitrogen content, total 
Chi (a + b) content, carotenoid content and greenness index from Fig. 3 



Effect 


Rubisco 


Soluble 


Total leaf 


Total 


Carotenoid 


Greenness 




content 


protein 


N 


Chi 




index 


Cultivar 


0.0018 


0.0079 


0.007 


<0.0001 


<0.0001 


<0.0001 


Growth CO2 


0.0127 


NS 


NS 


NS 


NS 


NS 


Stage 


<0.0001 


<0.0001 


<0.0001 


<0.0001 


<0.0001 


<0.0001 


Cultivar x CO2 


0.0361 


0.0048 


0.0026 


0.0289 


0.0022 


NS 


Cultivar x Stage 


0.0004 


0.002 


<0.0001 


0.0367 


0.0577 


NS 


CO2 X Stage 


NS 


NS 


NS 


NS 


0.0306 


NS 


Cultivar x CO2 x Stage 


NS 


NS 


NS 


NS 


NS 


NS 



Data are from the 2012 growing season. 

P-values <0.05 are in bold; P-values >0.05 but <0.10 are in italics; P-values >0.10 are indicated by NS. 



higher in Takanari than in Koshihikari at mid GF due to a decrease 
in gm in the latter. Flag leaves in the FACE treatment generally 
showed lower grn than in the ambient treatment in both years. 

Stomatal limitation to photosynthesis (SL) and mesophyll 
conductance limitation to photosynthesis (MCL) were calcu- 
lated for each CO2 response curve at the growth [CO2] 
(Ca = 390 and 590 |imol moP^ for ambient and FACE-grown 
plants, respectively) using Equations 1 and 2 (see the Materials 
and Methods and Supplementary Fig. S2 for details). Values of 
SL ranged from 4.2% to 8.6% over the mid and late GF stages in 
2012, and from 5.0% to 13.1% over the full heading and mid GF 
stages in 2013 (Fig. 2D). There was a significant overall cultivar 
effect for SL in both years (Table 1 ). SL was significantly lower in 
Takanari than in Koshihikari in the ambient treatment at late 
GF in 2012 and full heading and mid GF in 2013 (Fig. 2D). SL in 



the FACE treatment was significantly lower than ambient treat- 
ment in both years (Fig. 2D). 

Values of MCL ranged from 12.7% to 30.6% in 2012 and from 
1 5.6% to 26.7% in 201 3. Overall, MCL was lower in FACE than in the 
ambient treatments (Fig. 2E). There was no difference in MCL 
between cultivars at mid GF in 2012 and full heading in 2013, 
but MCL was significantly lower in Takanari than in Koshihikari 
at late GF in 2012 for both CO2 treatments and at mid GF in 2013 
for the FACE treatment (Fig. 2E). 

Leaf protein, nitrogen and photosynthetic 
pigment content 

The uppermost fully expanded leaves were sampled at several 
time points from the booting stage (late July) to late GF stage 
(early September) in the 2012 growing season. For total Rubisco, 
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Fig. 2 Photosynthetic parameters derived from analysis of CO2 response curves of Koshihikari (white bars) and Takanari (gray bars) grown under 
ambient (plain bars) and FACE (cross-hatched bars) conditions. Data were taken in 2012 at the mid grain filling (GF) (2 weeks after heading) and 
late GF stage (4 weeks after heading), and in 2013 at full heading (HD) (3 d after heading) and the mid GF stage. All measurements were taken on 
flag leaves. Error bars are SEMs. Different letters above each bar indicate significant differences between means at a = 0.05 within each meas- 
urement time point. n = 4. (A) Maximum carboxylation rate of Rubisco (Vc^max)- (B) Maximum electron transport rate (Jmax)- (Q Mesophyll 
conductance to CO2 (gm)/ expressed in |imol m~^ s~^ Pa~^ when atmospheric pressure is 101.3 kPa, this is nearly equal to 10 times the value of 
gm expressed in mol m~^ s~\ (D) Stomatal limitation to photosynthesis (SL) at each treatment's respective growth [CO2]. (E) Mesophyll 
conductance limitation to photosynthesis (MCL) at each treatment's respective growth [CO2]. 



leaf soluble protein and N content, there was a clear effect of 
growth stage on the values in both cultivars (Fig. 3A-C; Table 2). 
Values increased and peaked shortly after heading and steadily 
decreased afterwards. In both cultivars, the peak values of 
Rubisco, soluble protein and leaf N were observed at early GF 
(1 week after heading), but Takanari had significantly higher 
values than Koshihikari at mid GF (2 weeks after heading) 
(Fig. 3A-C). 

There was a significant cultivar by growth CO2 interaction 
effect for Rubisco, soluble protein and total leaf N content 
across all the time points as a whole (P = 0.036, 0.005 and 
0.003, respectively; Table 2). In Koshihikari, Rubisco, soluble 
protein and leaf N content were significantly decreased under 
the FACE treatment over the entire sampling period. In con- 
trast, in Takanari, there was no difference in Rubisco, total 
soluble protein and leaf N content between ambient- and 
FACE-grown leaves. 

Total Chi (Chi c? + b) and carotenoid content was signifi- 
cantly higher in Takanari than in Koshihikari across the entire 
sampling period (Fig. 3D, E; Table 2). On average, Takanari had 
29.3% and 45.5% higher Chi content than Koshihikari in the 
ambient and FACE treatments, respectively. In both cultivars, 
the Chi and carotenoid content reached a peak at mid GF and 
decreased thereafter. 



The greenness index, an indicator of leaf senescence, was 
about 4.3-4.6 for both cultivars from booting stage until mid 
GF, after which it fell sharply (Fig. 3F). From the mid-late to late 
GF stage (3 and 4 weeks after heading), the greenness index was 
higher in Takanari than in Koshihikari. There was no significant 
effect of growth [CO2] on greenness index (Table 2). 

The relationships of Rubisco, leaf soluble protein and Chi to 
total leaf N were plotted in Fig. 4. All the relationships showed a 
significant positive linear correlation, and plants grown under 
both ambient and FACE conditions fell along the same line for 
each cultivar. The relationships of Rubisco content vs. leaf N 
and leaf soluble protein vs. leaf N were essentially the same for 
both cultivars (Fig. 4A). However, Takanari showed a signifi- 
cantly higher Chi content than Koshihikari at the same leaf N 
level (Fig. 4B; P = 0.021). 



Discussion 




The results indicate that the source capacity in the flag leaves of 
Takanari remains higher than that of Koshihikari even when 
grown in elevated [CO2] in the field. Under ambient conditions, 
the overall leaf carbon assimilation rate, stomatal conductance, 
biochemical capacity of photosynthesis (V^max and jmaxX Chi 
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Growth stage 



Growth stage 



Growth stage 



Fig. 3 Leaf protein, nitrogen and pigment content of Koshihikari (white bars) and Takanari (gray bars) grown under ambient (plain bars) and 
FACE (cross-hatched bars) conditions. Flag leaves were sampled in 2012 at the following five growth stages: (1) booting (1.5 weeks before 
heading); (2) early grain filling (GF) stage (1 week after heading); (3) mid GF stage (2 weeks after heading); (4) mid-late GF stage (3 weeks after 
heading); and (5) late GF stage (4 weeks after heading). Error bars are SEMs. N = 4. (A) Rubisco content. (B) Total leaf soluble protein. (C) Total 
leaf nitrogen. (D) Total Chi (a + b). (E) Carotenoid content. (F) Greenness index (total Chl/carotenoid). 



content and carotenoid content were higher in Takanari com- 
pared with Koshihikari on a leaf area basis (Figs. 1-3). The 
varietal difference in stomatal conductance was particularly 
striking, as the gs of Takanari was at times as high as approxi- 
mately double that of Koshihikari (Fig. IB). These results are 
consistent with earlier reports (Hirasawa et al. 2010, Taylaran 
et al. 2011), and reflect the high-yield, high-productivity pedi- 
gree of Takanari (Takai et al. 2012, Takai et al. 2013). This ad- 
vantage of Takanari over Koshihikari was maintained when the 
plants were grown in elevated [CO2] under FACE conditions in 
both 2012 and 2013. Growth in eC02 caused an enhancement 
in the net carbon assimilation rate of both cultivars to a similar 
degree (Fig. lA). 

Higher stomatal conductance in Takanari over Koshihikari 
in ambient conditions resulted in significantly lower overall SL 
than in Koshihikari in both years (Table 1). The absolute value 
of SL in the ambient CO2 treatment ranged from 6.5% to 13.1% 
across both years, but the cultivar difference in SL at any given 
growth stage was only at most 2.2% and 4.0% in the ambient 
CO2 treatment in 2012 and 2013, respectively (Fig. 2D). In the 
FACE treatment, the absolute value of SL was lower than am- 
bient, ranging from 4.2% to 10.1%, and the maximum cultivar 
difference of SL in FACE was only 1.6% and 2.3% in 2012 and 
2013, respectively. 



In contrast, the limitation to photosynthesis due to CO2 
diffusion from the intercellular airspace to the site of carboxyl- 
ation (i.e. the MCL) was relatively higher than the SL, ranging 
from 12.7% to 30.6 % (Fig. 2E). The limitation of photosynthesis 
by mesophyll conductance was highest at the late GF stage, and 
it is also here that we observed the biggest cultivar difference in 
MCL; Takanari had about 36% higher across both ambient 
and FACE treatments, resulting in about 7.5% lower MCL than 
Koshihikari in both ambient and elevated [CO2] growth 
conditions (Fig. 2E). The importance of g^ in constraining 
photosynthesis has been explored extensively in the literature 
(e.g. Niinemets et al. 2009, Adachi et al. 2013), and the ability of 
Takanari to maintain higher g^ than Koshihikari late in 
the season may enable relatively higher rates of leaf photosyn- 
thesis even when photosynthetic capacity (V^max and JmaxX 
Rubisco and total leaf N are similar between cultivars 
(Figs. 2A, B, 3A, C). 

In comparison with the ambient [CO2] treatment, growth 
under FACE conditions resulted in lower values of total leaf N, 
leaf soluble protein content and Rubisco content in Koshihikari 
as early as the booting stage, and this difference increased as 
development progressed (Fig. 3A-C). In contrast, Takanari 
showed no significant difference in leaf N between ambient 
and elevated [CO2] treatments and consequently almost no 
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Fig. 4 Relationship between total leaf N and Rubisco N, soluble pro- 
tein N or Chi content, constructed from the time course data in Fig. 3. 
N content is assumed to be 16% of total protein. (A) Scatterplots of 
Rubisco N content vs. total leaf N and leaf soluble protein N vs. total 
leaf N. Circles indicate total soluble protein data, and triangles indicate 
Rubisco content data. Koshihikari is represented by open symbols, and 
Takanari is represented by filled symbols. (B) Scatterplot of total Chi 
content vs. total leaf N. Regression lines in (A) are drawn using pooled 
data from both cultivars for both ambient- and elevated [CO2] -grown 
treatments (n = 20); in (B), separate regression lines are drawn for each 
cultivar (n - 10). Error bars are SEMs of each treatment. 

change in leaf soluble protein content. Since leaf N content is 
one of the primary determinants of leaf photosynthesis (Evans 
1989, Makino 2003), this could be one of the factors driving the 
ability of Takanari to achieve very high photosynthetic rates 
while growing under FACE conditions. It is possible that leaf N 
content does not decrease in Takanari when the plants are 
grown in eC02 because of its heavy investment in root 
growth, especially from the panicle initiation stage onwards, 



and thus possesses a greater capacity for N uptake than 
Koshihikari during heading (Taylaran et al. 2011). However, 
this needs to be confirmed through future investigation of 
the total N uptake and N allocation within the whole plant. 

The largest difference between cultivars in Rubisco, soluble 
protein and total leaf N was observed at the mid GF stage 
(2 weeks after heading) (Fig. 3A-C), and this coincided with 
higher Vcmax ^nd J^ax and higher midday photosynthetic rates 
in Takanari than in Koshihikari at that growth stage in both 
2012 and 2013 (Figs. 1A, 2A, B). This suggests that Takanari 
may be able to maintain a higher leaf N status and thus higher 
photosynthetic capacity further into its GF period than 
Koshihikari, which would contribute to greater overall effi- 
ciency of GF before maturation (Nagata et al. 2001). In fact, 
Takanari shows the ability to fill its inferior spikelets just as well 
as its superior spikelets, in contrast to Koshihikari (G. Zhang, 
unpublished data), and our results here indicate a potential 
mechanism of longer duration of carbon assimilate supply 
during GF to enable that phenomenon. By the mid-late and 
late GF stages (3 and 4 weeks after heading), the leaf N, soluble 
protein and Rubisco content of both ambient-grown and 
FACE-grown Takanari were similar to those of FACE-grown 
Koshihikari (Fig. 3A-C). Ambient-grown Koshihikari showed 
higher values at the last sampling time point, perhaps indicating 
a slowing down of nutrient translocation into the grain and 
presence of sink capacity limitation. 

Takanari showed a significantly higher pigment content (Chi 
a + b and carotenoids) than Koshihikari from as early as the 
booting stage until near the end of GF (Fig. 3D, E). The Chi 
content vs. leaf N relationship showed a consistently higher Chi 
content in Takanari over Koshihikari at the same leaf N content 
(Fig. 4B), indicating that Takanari invests more of its N in light 
capture, although the impact of this difference on actual photo- 
synthetic rates is not clear. 

Why does Takanari keep its flag leaves greener than 
Koshihikari, even until the late stages of GF (Fig. 3D)? One 
explanation might be found in the canopy architecture of 
Takanari, which is distinctly different from that of Koshihikari 
(Xu et al. 1997b). In Koshihikari, the orientation of the flag 
leaves is relatively horizontal, and the panicles remain either 
above or at the same height as the flag leaves after heading, 
resulting in shading of the flag leaves and lower total intercep- 
tion of radiation as GF progresses (Fig. 5, left side). However, in 
Takanari, the flag leaves are comparatively wider and taller, and 
maintain a vertical orientation throughout the GF stage (Fig. 5, 
right side). This ensures that its flag leaf continues to be 
exposed to full light even as its large panicles become heavy 
and sag downwards during GF, and the presence of Chi may 
enable continued photosynthesis and carbohydrate production 
even at the later stages of GF, resulting in a high efficiency of GF 
(Xu et al. 1997a, Nagata et al. 2001, Saitoh et al. 2002). 

There is very limited intraspecific variability in the properties 
of Rubisco in rice, so the potential to increase rice productivity 
through breeding for direct improvements in the leaf biochem- 
ical properties may be small (Makino 2011). Instead, increased 
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Fig. 5 A colour photograph showing the canopies of Koshihikari (left) 
and Takanari (right) during the mid-grain filling stage. 



stomatal and/or mesophyll conductance, more efficient canopy 
structure or better optimized within-piant N allocation may 
be more useful targets for breeding higher productivity in 
rice (Saitoh et al. 2002, Hirasawa et al. 2010, Makino 2011, 
Adachi et al. 2013). Our results indicate that Takanari per- 
formed well in all of these aspects when grown in elevated 
[CO2] under field conditions, and would be a useful genetic 
resource for breeding new varieties of higher yielding rice that 
can thrive in future as well as present atmospheric conditions. 



until September 13; in 2013, fumigation ran from May 28 to 
September 13. In both years, 50% heading for Koshihikari was 
achieved by August 3, and for Takanari by August 9 and 6 in 
2012 and 2013, respectively. 

Gas exchange measurements 

Measurements of photosynthesis and stomatal conductance 
were conducted using a portable photosynthetic gas exchange 
system (LI-6400, Licor Inc.). Measurements were taken on the 
uppermost fully expanded leaves (two leaves for each cultivar 
per experimental plot per time point) at a photosynthetic 
photon flux density of 1,500-1,800 |imol m~^ s~^ and leaf 
chamber temperature of 32°C. The average leaf temperature 
within the chamber was 34.1±0.9°C and 33.1±0.5°C for 
Koshihikari and Takanari, respectively. The ambient air in the 
rice paddy was typically very humid, so the humidity was con- 
trolled through the occasional use of desiccant (Drierite, W.A. 
Hammond Drierite Co. Ltd.), resulting in a relative humidity 
within the leaf chamber of 79.4 ±5.6% for Koshihikari and 
80.5 ±2.4% for Takanari. Vapor pressure deficit was 
1.36 ± 0.34 kPa for Koshihikari and 1.08 ± 0.14 kPa for Takanari. 

CO2 response curves, consisting of 12 levels of CO2 ranging 
from 50 to 1,500 |imol mor\ were measured on two leaves per 
cultivar per ring during the mid GF and late GF growth stages in 
2012 (approximately 2 and 4 weeks after each cultivar's respect- 
ive heading dates), and during the full heading and mid GF 
stage in 2013 (3 days and 2 weeks after heading, respectively). 



Materials and Methods 
Site description 

We conducted the study at the Tsukuba FACE experimental fa- 
cility in Tsukubamirai City, Ibaraki Prefecture, Japan (35°58^N, 
139°60^E, 10 m a.s.l.) in the summers of 2012 and 2013. The soil 
is a Fluvisol, which is typical of alluvial areas. The average air tem- 
perature at the site in 2012 was 19.8°C in June, 24.9°C in July, 
26.8°C in August and 25.5°C for the first half of September; in 
2013, it was 21.3, 25.0, 26.9 and 25.1°C, respectively. The average 
ambient [CO2] at the site across the entire growing season (June- 
September) was 383.1 |imol moP^ in 2012 and 383.7 |imol moP^ 
in 2013. There were four control (ambient) plots and four FACE 
plots at the site (n = 4). The target concentration of the elevated 
[CO2] treatment (hereafter referred to as FACE treatment) was 
200 |imol moP^ above ambient [CO2], and the actual season-long 
mean [CO2] in the FACE plots was 577.1 |imol moP^ in 2012 and 
576.0 |imol moP^ in 2013. Further details regarding the experi- 
mental site set-up and CO2 control performance can be found 
in Nakamura et al. (2012), and soil chemical properties are given in 
detail in Hasegawa et al. (2013). 

Two rice {Oryza saiwa L.) cultivars were used in this study. 
'Koshihikari' is a japonica variety, and 'Takanari' is an indica 
variety. Three-week-old seedlings were transplanted into the 
experimental plots on May 23-24 in 2012 and May 22-23 in 
2013. CO2 fumigation was begun in 2012 on May 30 and ran 



Estimation of photosynthetic parameters and 
mesophyll conductance 

According to the Farquhar et al. (1980) steady-state biochemical 
model of leaf photosynthesis, the biochemical capacity for 
photosynthesis can be defined by two key parameters, namely, 
^cmax ^rid Jrnax (Farquhar et al. 1980, von Caem merer 2013). 
A third biochemical limitation to photosynthesis may also 
occur in the form of triose-phosphate utilization limitation 
(Sharkey 1985), but will not be considered here due to the in- 
ability to confirm conclusively its presence in the field study. 
These parameters can be estimated via curve fitting of the rela- 
tionship between the net carbon assimilation rate and the mole 
fraction of CO2 within the intercellular airspace of the leaf (Cj) 
[detailed equations may be found in Long and Bernacchi (2003) 
and von Caemmerer (2013)]. An Excel Solver-based non-linear 
curve-fitting utility based on these equations (Sharkey et al. 2007) 
was used simultaneously to solve for Vcmax/ imax leaf 
temperature. This utility employs the curve-fitting method of 
Ethier and Livingston (2004) to estimate from gas exchange 
data. To eliminate temperature differences between individual 
leaves at the time of measurement, Vcmax and J^ax were adjusted 
to 32°C using the temperature response functions of Bernacchi 
et al. (2001, 2003) which have been implemented in the Sharkey 
et al. (2007) curve-fitting utility; values of mesophyll conductance 
were adjusted to 32°C using the temperature response function 
of Scafaro et al. (2011) for 0. sativa. 
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A-Cj data were also fit under the assumption that gm is 
infinite and Cj^Q to produce estimates of Vcmax and J^ax 
from the 'conventional' curve-fitting method (described by 
Long and Bernacchi 2003). We found that this method pro- 
duced a poorer fit of the measured data points, so only results 
from the curve-fitting method described above (Sharkey et al. 
2007) are used in this study. Example A-Q and A-Q fitted 
curves can be found in Supplementary Fig. SI, and the 
mean parameter estimates derived from the conventional 
curve-fitting are also available (Supplementary Table SI). 

Calculation of CO2 diffusion limitations to 
photosynthesis 

Limitations to photosynthesis due to resistance in the CO2 
diffusion pathway from the atmosphere to the site of carboxyl- 
ation were calculated for each A-Q curve. This was divided into 
two components: limitation of photosynthesis due to CO2 re- 
sistance at the leaf boundary layer and stomata, SL; and limita- 
tion of photosynthesis due to resistance from the intercellular 
airspace to the site of carboxylation, MCL. 

SL (expressed as a percentage) was calculated following the 
method described in Long and Bernacchi (2003): 



SL = 



(Ao-Ai) 



X 100 



0) 



where Aq represents the theoretical rate of photosynthesis 
when we assume that the boundary layer conductance to 
CO2, stomatal conductance to CO2 and gm infinite, i.e. if 
Q = Q (atmospheric mole fraction of CO2, assumed to be 390 
and 590 jimol moP^ for the ambient and FACE-grown plants, 
respectively); Ai represents the theoretical rate of photosynthe- 
sis when only is assumed to be infinite, i.e. if Q = Q. 

MCL (expressed as a percentage) was calculated in a similar 
manner to SL, as follows: 



MCL : 



(At - A2) 



100 



(2) 



where is the same as above, and A2 represents the actual 
measured photosynthesis rate at when the measured sto- 
matal conductance to CO2 and estimated g^ are taken into 
account. A graphical representation of the calculation of SL and 
MCL in the context of the CO2 response curve can be found in 
Supplementary Fig. S2. 

Biochemical analyses 

The uppermost fully expanded leaves of both cultivars were 
sampled in 2012 at five time points from the booting stage 
(late July, approximately 1.5 weeks before heading) to the late 
GF stage (early September, 4 weeks after heading). Three leaves 
per cultivar were sampled in each experimental plot per time 
point. 

Approximately 6 cm^ of leaf material was sampled from the 
same leaves and immediately frozen in liquid nitrogen and 
stored at — 80°C. After thawing, 0.5 cm^ of leaf material was 



placed in 1.8 ml of 95% ethanol and stored in the dark at 4°C 
for 3d to extract the leaf pigments. The absorbance of the 
solutions at 470, 648.6 and 664.1 nm was then measured in a 
spectrophotometer (Model U-1800, Hitachi High-Technologies 
Corp.). The total Chi content, carotenoid content and green- 
ness index (the ratio of total Chi to carotenoid content) were 
calculated from these absorbance values using the equations of 
Lichtenthaler and Buschmann (2001). Another 1 cm^ of leaf 
material was ground using a glass mortar and pestle chilled at 
4°C with 500 |il of extraction buffer consisting of 50 mM Tris- 
HCI (pH 7.6), 5% glycerol, 2mM Na-iodoacetate, 5mM dithio- 
threitol and 0.2% (w/v) polyvinylpyrrolidone. Leaf soluble 
protein was quantified from a portion of this crude extract 
using a commercial version of the assay of Bradford (1976) 
(Quick Start Bradford Protein Assay, Bio-Rad Laboratories Inc.). 

A 30 |il aliquot of the crude protein extract was mixed with 
30 |il of 1% lithium dodecyl sulfate buffer and heated at 97°C for 
3min. Then 6|il of this mixture was loaded into a 7.5-15% 
polyacrylamide gradient gel and electrophoresis was conducted 
at 150 V for 30min using the XV Pantera MP electrophoresis 
system (DRC Co., Ltd.). The gel bands were stained using Quick- 
CBB Plus (Wako Pure Chemical Industries, Ltd.). 

Densitometry analysis was performed to estimate the quan- 
tity of Rubisco large subunit. Total Rubisco content was calcu- 
lated as 1.24 times the Rubisco large subunit content, based on 
the ratio of molecular weights of the large and small subunits. 
A single stock of bovine standard albumin (A7906, Lot # 
SLBC0647V, Sigma-Aldrich Co.) was used to create a set of 
protein standards to approximate the quantity of Rubisco pro- 
tein in each gel. 

Approximately 10 cm^ was taken from the middle section of 
the leaf blade and dried in an oven at 80°C for at least 3 d. The 
dried leaf tissue was used to quantify total leaf N content using 
an NC analyzer (Sumigraph NC-22, SCAS Ltd.). 

Statistical analysis 

Gas exchange and biochemical data were analyzed for each year 
as a randomized complete block split-plot design using a gen- 
eralized linear mixed model and three-way (cultivar x growth 
[CO2] X growth stage) analysis of variance (PROC GLIMMIX, 
SAS 9.3, SAS Inc.) to examine overall treatment effects, as well as 
a two-way (cultivar x growth [CO2]) analysis of variance at 
each growth stage to identify treatment mean differences. 
The FACE treatment was treated as a whole-plot effect, and 
the cultivars were subplots. 



Supplementary data 

Supplementary data are available at PCP online. 



Funding 

This work was supported by the Ministry of Agriculture, 
Forestry and Fisheries, Japan, [through a research project 



PteCe/ZPfiys/o/. 55(2): 381-391 (2014) doi:10.1093/pcp/pcu009 © The Author 2014. 389 



J-^ ^ J-^ C p. Chen et al. 

PLANT & CELL PHYSIOLOGY 



entitled 'Development of technologies for mitigation and 
adaptation to climate change in Agriculture, Forestry and 
Fisheries']; the Japan Society for the Promotion of Science [a 
Grant'in-Aid for Scientific Research on Innovative Areas 
(No.24114711) as part of the project entitled 'Comprehensive 
studies of plant responses to high CO2 world by an innovative 
consortium of ecologists and molecular biologists']. 



Acknowledgments 




We thank the Hayatsu research group at NIAES for the gener- 
ous use of their laboratory facilities for the biochemical analyses. 
We also thank Professor Amane Makino of Tohoku University 
for advice regarding the extraction of leaf protein and estima- 
tion of Rubisco content, and Guoyou Zhang for his assistance 
with leaf sampling. 



Disclosures 




The authors have no conflicts of interest to declare. 



References 




Adachi, S., Nakae, T., Uchida, M., Soda, K., Takai, T., Oi, T. et al. (2013) 
The mesophyll anatomy enhancing C02 diffusion is a key trait for 
improving rice photosynthesis, j. Exp. Bot 64: 1061-1072. 

Ainsworth, E.A. and Long, S.P. (2005) What have we learned from 15 
years of free-air C02 enrichment (FACE)? A meta-analytic review of 
the responses of photosynthesis, canopy. New Phytol. 165: 351-371. 

Ainsworth, E.A., Rogers, A. and Leakey, A.D.B. (2008) Targets for crop 
biotechnology in a future high-C02 and high-O-3 world. Plant 
Physiol. 147: 13-19. 

Alexandratos, N. and Bruinsma, J. (2012) World Agriculture Towards 
2030/2050: The 2012 Revision FAO, Rome. 

Bernacchi, C.J., Pimentel, C. and Long S.P. (2003) In vivo temperature 
response functions of parameters required to model RuBP-limited 
photosynthesis. Plant Cell Environ. 26: 1419-1430. 

Bernacchi, C.J., Singsaas, E.L, Pimentel, C, Portis, A.R. and Long, S.P. 
(2001) Improved temperature response functions for models of 
Rubisco-limited photosynthesis. Plant Cell Environ. 24: 253-259. 

Bradford, M.M. (1976) A rapid and sensitive method for the quanti- 
tation of microgram quantities of protein utilizing the principle of 
protein-dye binding. Anal. Biochem. 72: 248-254. 

Ethier, G.J. and Livingston, N.J. (2004) On the need to incorporate 
sensitivity to C02 transfer conductance into the Farquhar-von 
Caem merer- Berry leaf photosynthesis model. Plant Cell Environ. 
27: 137-153. 

Evans, J.R. (1989) Photosynthesis and nitrogen relationships in leaves 
of C3 plants. Oecologia 78: 9-19. 

Farquhar, CD., Caemmerer, S.V. and Berry, J.A. (1980) A biochemical- 
model of photosynthetic CO2 assimilation in leaves of C-3 species. 
Planta 149: 78-90. 

Hasegawa, T, Sakai, H., Tokida, T., Nakamura, H., Zhu, C, Usui, Y. et al. 
(2013) Rice cultivar responses to elevated C02 at two free-air 
C02 enrichment (FACE) sites in Japan. Funct. Plant Biol. 40: 
148-159. 



Hatfield, J.L, Boote, K.J., Kimball, B.A., Ziska, L.H., Izaurralde, R.C., 
Ort, D. et al. (2011) Climate impacts on agriculture: implications 
for crop production. Agron. j. 103: 351-370. 

Hendrey, G.R. and Miglietta, F. (2006) FACE technology: past, present, 
and future. In Managed Ecosystems and C02: Case Studies, 
Processes, and Perspectives. Ecological Studies Vol. 187. Edited by 
Nosberger, J., Long, S.P., Norby, R.J., Stitt, M., Hendrey, G.R. and 
Blom, H. pp. 15-43. Springer, Berlin. 

Hirasawa, T, Ozawa, S., Taylaran, R.D. and Ookawa, T. (2010) Varietal 
differences in photosynthetic rates in rice plants, with special 
reference to the nitrogen content of leaves. Plant Prod. Sci. 13: 
53-57. 

Imbe, T., Akama, Y., Nakane, A., Hata, T., Ise, K., Ando, I. et al. (2004) 
Development of a multipurpose high-yielding rice variety 'Takanari'. 
(In Japanese with English summary). Bull. Natl Inst. Crop Sci. 5: 35-51. 

Khush, G.S. (2005) What it will take to feed 5.0 billion rice consumers 
in 2030. Plant Mol. Biol. 59: 1-6. 

Leakey, A.D.B., Ainsworth, E.A., Bernacchi, C.J., Rogers, A., Long, S.P. 
and Ort, D.R. (2009) Elevated C02 effects on plant carbon, nitrogen, 
and water relations: six important lessons from FACE.j. Exp. Bot. 60: 
2859-2876. 

Lichtenthaler, H.K. and Buschmann, C. (2001) Chlorophylls and carot- 
enoids: measurement and characterization by UV-VIS spectroscopy. 
In Current Protocols in Food Analytical Chemistry (CPFA). Edited 
by Wrolstad, R.E., Acree, T.E., An, H., Decker, E.A, Penner, M.H., Reid, 
D.S. et al. pp. F4.3.1-F4.3.8. John Wiley & Sons, New York. 

Long, S.P. and Bernacchi, CJ. (2003) Gas exchange measurements, 
what can they tell us about the underlying limitations to photo- 
synthesis? Procedures and sources of error, j. Exp. Bot. 54: 
2393-2401. 

Makino, A. (2003) Rubisco and nitrogen relationships in rice: leaf 
photosynthesis and plant growth. So/7 Sci. Plant Nutr. 49: 319-327. 

Makino, A. (2011) Photosynthesis, grain yield, and nitrogen utilization 
in rice and wheat. Plant Physiol. 155: 125-129. 

Meehl, G.A., Stocker, T.F., Collins, W.D., Friedlingstein, P., Gaye, A.T, 
Gregory, J.M. et al. (2007) Global climate projections. In Climate 
Change 2007: The Physical Science Basis. Contribution of Working 
Group I to the Fourth Assessment Report of the Intergovernmental 
Panel on Climate Change. Edited by Solomon, S., Qin, D., Manning, 
M., Chen, Z., Marquis, M., Averyt, K.B. et al. Cambridge University 
Press, Cambridge. 

Nagata, K., Yoshinaga, S., Takanashi, J.-i. and Terao, T. (2001) Effects of 
dry matter production, translocation of nonstructural carbohy- 
drates and nitrogen application on grain filling in rice cultivar 
Takanari, a cultivar bearing a large number of spikelets. Plant 
Prod. Sci. 4: 173-183. 

Nakamura, H., Tokida, T., Yoshimoto, M., Sakai, H., Fukuoka, M. and 
Hasegawa, T. (2012) Performance of the enlarged Rice-FACE system 
using pure C02 installed in Tsukuba, Japan, j. Agric. Meteorol. 68: 
15-23. 

Niinemets, T, Diaz-Espejo, A., Flexas, J., Galmes, J. and Warren, C.R. 
(2009) Role of mesophyll diffusion conductance in constraining 
potential photosynthetic productivity in the field, j. Exp. Bot. 60: 
2249-2270. 

Nosberger, J., Long, S.P., Norby, R.J., Hendrey, G.R. and Blum, H. (2006) 
In Managed Ecosystems and CO2: Case Studies, Processes, and 
Perspectives Springer, Berlin. 

Saitoh, K., Yonetani, K., Murota, T. and Kuroda, T. (2002) Effects of flag 
leaves and panicles on light interception and canopy photosynthe- 
sis in high-yielding rice cultivars. Plant Prod. Sci. 5: 275-280. 



390 Plant Cell Physiol. SS{2): 38^-39^ (2014) doi:10.1093/pcp/pcu009 © The Author 2014. 



Takanari leaf physiology under FACE treatment ^-^^ ^ 

PLANT & CELL PHYSIOLOGY 



Scafaro, A.P., Von Caemmerer, S., Evans, J.R. and Atwell, B.J. (2011) 
Temperature response of mesophyll conductance in cultivated and 
wild Oryza species with contrasting mesophyll cell wall thickness. 
Plant Cell Environ. 34: 1999-2008. 

Sharkey, T.D. (1985) O-2-insensitive photosynthesis in C-3 plants— its 
occurrence and a possible explanation. Plant Physiol. 78: 71-75. 

Sharkey, T.D., Bernacchi, C.J., Farquhar, CD. and Singsaas, E.L (2007) 
Fitting photosynthetic carbon dioxide response curves for C-3 
leaves. Plant Cell Environ. 30: 1035-1040. 

Takai, T., Adachi, S., Taguchi-Shiobara, F., Sanoh-Arai, Y., Iwasawa, N. 
et al. (2013) A natural variant of NAL1, selected in high-yield rice 
breeding programs, pleiotropically increases photosynthesis rate. 
Sci. Rep. 3: 2149. 

Takai, T., Arai-Sanoh, Y., Iwasawa, N., Hayashi, T., Yoshinaga, S. and 
Kondo, M. (2012) Comparative mapping suggests repeated 
selection of the same quantitative trait locus for high leaf photo- 
synthesis rate in rice high-yield breeding programs. Crop Sci. 52: 
2649-2658. 

Tausz, M., Tausz-Posch, S., Norton, R.M., Fitzgerald, G.J., Nicolas, M.E. 
and Seneweera, S. (2013) Understanding crop physiology to select 
breeding targets and improve crop management under increasing 
atmospheric C02 concentrations. Environ. Exp. Bot. 88: 71-80. 

Taylaran, R.D., Adachi, S., Ookawa, T., Usuda, H. and Hirasawa, T. 
(2011) Hydraulic conductance as well as nitrogen accumulation 



plays a role in the higher rate of leaf photosynthesis of the most 
productive variety of rice in Japan, j. Exp. Bot. 62: 4067-4077. 
Taylaran, R.D., Ozawa, S., Miyamoto, N., Ookawa, T., Motobayashi, T. 
and Hirasawa, T. (2009) Performance of a high-yielding modern rice 
cultivar Takanari and several old and new cultivars grown with and 
without chemical fertilizer in a submerged paddy field. Plant Prod. 
Sci. 12: 365-380. 

von Caemmerer, S. (2013) Steady-state models of photosynthesis. 
Plant Cell Environ. 36: 1617-1630. 

Xu, Y.F., Ookawa, T. and Ishihara, K. (1997a) Analysis of the photo- 
synthetic characteristics of the high-yielding rice cultivar Takanari. 
jpn. j. Crop Sci. 66: 616-623. 

Xu, Y.F., Ookawa, T. and Isihara, K. (1997b) Analysis of the dry matter 
production process and yield formation of the high-yielding rice 
cultivar Takanari from 1991 to 1994. jpn. j. Crop Sci. 66: 42-50. 

Yadav, S.S., Redden, R., Hatfield, J.L, Lotze-Campen, H. and Hall, A.J.W. 
(2011) Crop Adaptation to Climate Change. Wiley. 

Ziska, LH., Bunce, J.A., Shimono, H., Gealy, D.R., Baker, J.T., 
Newton, P.C.D. et al. (2012) Food security and climate change: on 
the potential to adapt global crop production by active selection to 
rising atmospheric carbon dioxide. Proc. R. Soc. B: Biol. Sci. 279: 
4097-4105. 



PteCe/ZPfiys/o/. 55(2): 381-391 (2014) doi:10.1093/pcp/pcu009 © The Author 2014. 391 



